Shell evolution and its indication on the isospin dependence of the spin-orbit splitting 
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The available experimental data on shell evolution indicate that the strength of the spin-orbit 
(SO) single-particle potential may be enhanced in neutron-rich nuclei. We observe that such a 
simple scheme destroys the Harmonic Oscillator (HO) magic numbers N = 8 and 20 and generates 
new SO magic numbers like N = 6, 14, 16, 32 and 34. The traditional SO magic numbers like 
N = 28 and 50 and N = 14 seen in 22 O are eroded somehow in neutron-rich nuclei due to the 
sensitivity of larger-/ orbitals to the depth of the central potential but they are more robust than 
the HO magic numbers. The N = 82 shell closure persists in neutron-rich nuclei while HO shell 
closures like N — 40 and 70 do not emerge. Both mechanisms contribute to enhancing the N = 56 
and 90 gaps by splitting the lds/2 and 0(^/2 and the 0hg/ 2 and I/7/2 orbitals. 
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One of the most important and challenging frontiers of 
nuclear structure physics is the study of nuclei at the limit 
of stability, especially neutron-rich nuclei with weakly 
bound neutrons. A topic of particular interest is the evo- 
lution of the shell structure in those nuclei. That is, the 
magic number may change dramatically depending on 
the NJZ ratio when we move towards the particle drip 
lines Such study is important not only due to the 
expected variation in properties of nuclei and the forma- 
tion of island of inversion but also for the understanding 
of nuclear astrophysics as well as the nucleon-nucleon in- 
teraction. Nowadays it is rather commonly accepted that 
the N = 8 and 20 Harmonic Oscillator (HO) shell clo- 
sures disappear in neutron-rich nuclei [l|. On the other 
hand, new magic numbers like N = 14, 16, and 32 may 
emerge. Shell model calculations suggest that N = 34 
may also be a magic number in Ca isotopes, depending 
on the effective interactions used (see, e.g., Fig. 2 in Rcf. 
1). 

Among the intensive discussions on the mechanisms 
behind the shell evolution phenomena, we mention the 
tensor component of the two-body shell model interac- 
tion and the three-body interaction (see Refs. 0, Q and 
references therein). In that scheme, the evolution of the 
shells are solely determined by the correlation between 
valence nucleons in the open shell. The shell evolution 
has also been analyzed from a self-consistent mean-field 
point of view without assuming an inert core (see, e.g., 
Ref. 0, [![). It suggests that a significant part of the 
spin-orbit splitting may come from the two-body spin- 
orbit (SO) and tensor forces and three-body forces. The 
availability of experimental data in nuclei with large N/Z 
ratios may provide a ground to constrain the properties 
of different components of the interaction, such as the 
isovector channel of the SO interaction, which are not 
well defined but may be responsible for the shell evolu- 
tion. 

From a simple phenomenal point of view, the shell 



structure is characterized by the presence of gaps in 
single-particle spectra. The HO mean field with SO cou- 
pling was the first successful model to predict correctly 
the traditional magic numbers The calculated shell 
structure may change if an isospin dependence of the SO 
coupling is taken into account. To illustrate this point, 
in Fig. [T]we evaluate the HO single-particle spectra by 
adding an isospin-dependent SO coupling of the form 
A(l + K SO ; ^r^)^l-s H to the HO potential. A negative 
kso means that the SO coupling will gradually diminish 
when we go from N = Z to neutron-rich nuclei, which 
indicates that SO shell closures like N = 28 and N = 50 
will disappear and the spectra will be characterized by 
only HO magic numbers. However, this contradicts most 
available experimental information on the shell evolution 
even though there is indication that the N = 28 may 
have eroded in 42 Si 0413. On the other hand, a pos- 
itive Kso would indicate that the SO coupling will be 
enhanced in neutron-rich nuclei. In such a scheme, the 
A = 8 and 20 HO magic numbers will disappear while 
SO shell closures like N = 6, 14, 16, 32, 34 and 56 will 
appear. It should be kept in mind, however, that the 
gaps at N = 28 and 50 will be enhanced in neutron-rich 
nuclei within this naive picture. 

In reality, we may expect that the SO coupling will 
be reduced in neutron-rich drip-line nuclei with a dif- 
fusive surface since the SO interaction is peaked at the 
nuclear surface. One also has to consider that the single- 
particle orbitals may show different I dependence on the 
strength of the potential, which can lead to a systematic 
change of the shell structure [llj . Thus a more realistic 
description of the shell structure may be obtained with 
the Woods-Saxon potential which has a glorious history 
of success and is still one of the most suitable models in 
describing the nuclear single-particle structure. A variety 
of paramcterizations of the Woods-Saxon potential exists 
(see, e.g., Refs. gl [Tp-Hand Table II in Ref. [3). In 
the "standard" one [5], LL2| , the strengths of the central 
and SO potentials are given as 
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^ = ^o(l + ^t-T d ) 
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FIG. 1. (Color online) The evolution of the shell structure 
as a function of (N — Z)/A with the HO potential plus SO 
coupling of the form A(l + kso N ^ Z )fnu\ ■ s. We take A = 0.2 
and kso = — 1 (left) and 1 (right). The Ogg/2 orbital is shifted 
upwards by O.Shoj for a clearer presentation. 
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FIG. 2. (Color online) Evolution of the N=20, 28, 50 and 82 
gaps as a function of proton number Z for calculations with 
the standard Woods-Saxon parameter and kso — K (open 
circles) and — k (solid circles). The black squares denote the 
corresponding experimental data [l6| . 



and 

^o=AVo(l+^t-T d ), (2) 

where we have replaced the original N ~ Z term with 
4t ■ Td to get a consistent description of both protons and 
neutron orbitals. t and denote the isospin quantum 
numbers of the last nucleon and the daughter nucleus, 
respectively. The total isospin of the system is T = t + 
Tyt-i- We have 4t • T A _i = -3 for the T = ground 
state of a N = Z nucleus and 

4t • Ta-i = N — Z — 1 for neutron orbits 

= -(N - Z + 3) for proton orbits (3) 

in N > Z nuclei with T = {N - Z)/2 [T|. In Rcf. 
@, the isospin-dependent terms in Eqs. (1) and (2) are 
parameterized as 

k = nso = - — , (4) 

where the SO potential depth is assumed to have the 
same isospin dependence as that of the central potential. 
This assumption is rather commonly adapted (T3 . [Til ]. 
The typical strength of K is in the range —0.6 ~ —0.9. 

We have done a systematic calculations on the shell 
evolution with the standard Woods-Saxon parameters 
(see, also, Ref. [HI). It suggests that SO shell closures 
like N = 28, 50 and 82 will erode in neutron-rich nuclei 
but in a manner that may be too fast, as can be seen in 



Fig. [5] Moreover, it cannot explain in a straight-forward 
way the disappearance of HO shell closures like N = 8 
and 20. The predicted N = 40 gap is also too strong. 
As we will show below, it seems that the problem may 
be fixed if we assume a str ong po sitive kso- Calcula- 
tions with other parameters [l5[ will lead to the same 
conclusion. 

To get a qualitative idea on the role played by ksOi two 
kinds of calculations are done with the standard Woods- 
Saxon parameter Q by taking kso = k and —k. The 
numerical code GAMOW is used [T3|- Calculations on 
the evolution of the N=20, 28, 50 and 82 magic num- 
bers are plotted in Fig. [2j The TV = 20 shell closure 
is expected to disap pea r in neutron-rich nuclei like 32 Mg 
[l| (see, also, Ref. [HI). In calculations with the stan- 
dard parameter, however, the gap in nucleus like 28 is 
as large as 4.2 MeV. The N = 20 shell persist even if 
one takes kso = 0, as shown in Ref. [HJ. If one takes 
kso = —k, however, the N = 20 gap in 28 would reduce 
to only 1.2 MeV. 

The situation around N = 28 is somewhat compli- 
cated 0. The N — 28 will be reduced in neutron-rich 
nuclei as expected from the standard Woods-Saxon cal- 
culations But the reduction of this shell gap will 
be significantly retarded if a positive kso is assumed, 
as can be seen in Fig. [2] Similarly, the N = 82 shell 
closure may be reduced in very neutron-rich nuclei like 
110 Ni since the large I orbital 0/T. n / 2 will lose energy faster 
than other smaller-Z orbitals when the potential becomes 
shallower. Indeed, the distance between the 0/in/2 and 
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FIG. 3. (Color online) Evolution of the single-particle ener- 
gies of the IP3/2, lpi/2, O/5/2 and 0g 9 / 2 orbitals in Ca isotopes 
as a function of neutron number iV for calculations with the 
standard Woods-Saxon parameter and kso = k (open sym- 
bols) and — k (solid symbols). In former calculations, the 
N = 40 shell gap is significantly enhanced in neutron-rich 
isotopes while the N = 34 sub-shell gap is eroded. 



I/7/2 orbitals decrease from 4.1 MeV in 132 Sn to 0.65 
MeV in 110 Ni in calculations with the standard parame- 
ter. The gap is increased to 2.5 MeV if we simply take 
kso — ~ K - It should be mentioned that experimental 
information in Ni isotopes is known only up to V = 50. 
The I/7/2 orbital in 110 Ni is predicted to be unbound in 
both calculations. 

The sign of kso can have significant influence on the 
relative strength of the N = 40 and 50 shell gaps in 
neutron-rich nuclei. As examples, in Fig. [3] we plot the 
single-particle spectra of Ca isotopes calculated with the 
standard parameter. The N — 40 shell cap in 60 Ca is 
predicted to be as large as 4.4 MeV in calculations with 
the standard parameter. While in the unbound nucleus 
70 Ca the N = 50 gap is only 0.67 MeV. The N = 50 
shell gap may disappear if the Ogg/2 orbital lose too much 



energy due to the shallowing of the central potential [11 



On the other hand, the gaps are calculated to be 0.43 
MeV and 2.3 MeV, respectively, if we take kso = — k, 
which means that a positive kso will restore the N = 50 
magic number by "destroying" the N = 40 HO shell 
closure. There is no experimental indication that TV = 40 
will emerge as a new magic number in neutron-rich nuclei 

MM- 

The N = 32 gap between orbitals p 3 / 2 and px/2 is n °t 
much affected by the kso term. The N = 32 gap in 52 Ca 
is calculated to be 1.8 MeV with the standard parameter. 
It is increased to 2.3 MeV if we take kso = —k. The 
difference between lpi/ 2 and 0/ 5 / 2 is only 0.92 MeV in 
54 Ca. The N = 34 gap increased to 1.7 MeV if we take 
kso = -k. 



Recent experiments suggest that 22 and 24 should 
be doubly magic nuclei ]l| (see, also, Refs. [2lT[23l | for 
recent results). 24 O is the heaviest bound oxygen isotope 
that has been observed so far. The N = 16 gap between 
the neutron ls±/2 and states are measured to be 

4.86 ±0.13 MeV in Ref. 
lations with kso = k and 



21] . The values given by calcu- 
k are 3.3 MeV and 4.5 MeV, 
respectively. The calculated N = 14 gap in 20 O increases 
from 1.8 MeV to 2.7 MeV if one takes k S o = -«• 

The N = 14 gap disappear in C and N isotopes with 
nearly degenerate and IS1/2 orbitals. This is eas- 

ily understood since the Ods/2 orbital loses its energy 
faster when going towards the dripline, resulting a nearly- 
degenerate 0d 5 /2 and lsi/2 ll|. From a shell model point 



of view, where the above-mentioned mechanism is miss- 
ing, this fact is related to the complicated interplay be- 
tween the isovector and isoscalar two-body interactions 

M- 

A possible different form of isospin dependence in the 
SO potential than that of the central pot ential has been 
the subject of several studies [H, li~5l. |25|. Ref. [l5| as- 
sumed that kso — 0. The study of Ref. [l~3| showed 
that kso ~ 0.2 to 0.7, which is opposite to K, can also 
explain the single-particle spectra in neutron-rich nuclei 
132 Sn and 208 Pb. They suggest that such a opposite value 
is consistent with the two-body SO interaction as well as 
the Walecka and Skyrme-Hartree-Fock calculations. The 
study of Ref. showed that the isoscalar SO coupling in 
the Skyrme force is reduced while the tensor coupling is 
strongly attractive, which may also indicate that the SO 
splitting can be enhanced in neutron-rich nuclei. How- 
ever, it should be mentioned that the spectra of heavy 
stable nuclei are quite insensitive to the sign of kso- As 
a result, it_is_not determined by a normal global fitting 
procedure 



25] 



The effect of the kso term may be partly swallowed by 
other Woods-Saxon parameters, in particular the SO ra- 
dius parameter rso m the Woods-Saxon potential, rso 
was taken as a free parameter in several calculations 
[Til . (HI, which tends to take values smaller than that 
of the central potential, i.e., r$. To explore this point 
further, we re-fitted the Woods-Saxon parameters under 
the restriction kso = ~ K - The parameters of the Woods- 
Saxon potential are adjusted to single-particle and single- 
hole states around the doubly-magic nuclei 16 , 40,48 Ca, 
56 Ni, 100 Sn, 132 Sn and 208 Pb, as listed in Refs. MM- 
We use the same fitting procedure as employed in Ref. 
[26| . The fitting to nuclear single-particle energies usually 
favors larger values for the radius parameter rrj, which 
may lead to a bad description of nuclear charge radii and 
moments of inertia [27|. In this work we restrict that 
ro < 1.3 fm. The results thus obtained are presented in 
Table [T] Calculations on the evolution of the A" = 20, 28, 
50 and 82 magic numbers arc plotted in Fig. 2J which 
show similar trends to those in Fig. O 

In Fig. [5] we evaluate the influence of the isospin de- 
pendence in the SO coupling on the proton single-particle 
spectra in neutron-rich nuclei. As seen from Eq. (3), a 



4 



TABLE I. Woods-Saxon potential parameters obtained by fitting to single-particle and single-hole states around doubly-magic 
nuclei with the restriction kso = — k and comparison with some existing parameters. 
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0.647 


KSO = 


K 


Ref. [14] 


49.6 


1.347(n)/1.275(p) 1.31(n)/1.32(p) 
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FIG. 4. (Color online) Same as Fig. 
with the re-fitted parameter and kso 
kso = — k (solid symbols). 



[2] but for calculations 
= k (open symbols) and 
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FIG. 5. (Color online) Same as Fig. [4] but f° r calculations on 
the Z = 8, 20, 28 and 50 proton shell gaps as a function of 
neutron number N in neutron-rich nuclei. The open and solid 
symbols stand for calculations with kso = K and kso = —n, 
respectively. 



positive kso would suggest that the proton SO splitting 
is reduced in neutron-rich nuclei. Ref. (28| does show 
that the splitting between the proton shells Ogj/2 and 
increases with neutron excess in neutron- rich Sb 
isotopes. As already shown in Ref. [29| , this fact can be 
reproduced if we take a positive value for kso ■ Whereas 
it is related the effect of the two-body tensor force within 
the Skyrme-Hartree-Fock approach [30j- It is difficult 
to conclude on the situation in light neutron-rich nuclei 
since experimental results are still inadequate. However, 
shell-model calculations tend to suggest that the split- 
tings between SO partners like 0p 3 / 2 and Op-, /? and O/7/2 
and O/5/2 may reduce with neutron excess 

gHU. This 

is consistent with the binding energy systematics in Ref. 
[l[ which shows that the Z = 8 and 20 gaps increase in 
neutron-rich nuclei. On the other hand, the Z = 28 gap 
may decrease 

[Elm. A11 these facts 

are consistent with 

our assumption that kso is positive. 

In this work we concentrated on the structure of light 
nuclei. But as shown in Ref. [l|j], the choice of a posi- 
tive Kso is consistent with the single-particle spectra in 
132 Sn and 208 Pb. Indeed, calculations with the re-fitted 
parameter in Table Q] can describe these states equally 
well in c omp arison with those with the potentials given 
in Refs. [M0J]. 

The single-particle spectra may be influenced by cor- 
relation effects including deformation, particle vibration 
and pairing [27j . To explore the influence of pairing cor- 
relation, the single-particle plus pairing Hamiltonian is 



FIG. 7. (Color online) Potential energy surface (PES) calcu- 
lations for 32 Mg ground state with the re-fitted Woods-Saxon 
paramter and kso = — K (left) and k (right). For each (,82,7) 
point the energy is minimized with respect to the 04, deforma- 
tion. The interval between neighboring contours is 0.1 MeV. 



solved exactly with a Lanczos approach (33[ ■ We assume 
that the ground states of even-even nuclei are all paired 
with seniority zero whereas those of odd- A nuclei are as- 
sumed to be of seniority one. We take the semi-magic Ca 
isotopes as examples, which have been studied recently 
both experimental and theoretically 0, |34T - [36| . Calcu- 
lations are done within a model space containing the 
O/7/2, 0/ 5/2 , lp 3/2 , lpi/2, Ogg/2 and ld 5/2 neutron or- 
bitals by assuming 40 Ca as the core. To minimize the 
number of free parameters, a simple constant pairing 
strength, G = 1.795 MeV, is employed in all calcula- 
tions. The ground state energies thus calculated are pre- 
sented in Fig. [(^together with experimental data. Calcu- 
lations with kso = —n predict larger binding energies for 
neutron-rich nuclei 51_58 Ca than the extrapolation val- 



ues given in Ref. [37J. This is consistent with Ref. [34 1 
where additional binding is measured for nuclei 51 ' 52 Ca. 
The increased binding in Ca isotopes may indicate a sig- 
nicant subshell gap at TV = 32 [3J|. There is a kink in 
the systematics of calculated two-neutron separation en- 
ergies at TV = 34, which suggests that TV = 34 may also 
be a subshell. However, the lpi/2 particle and IP3/2 hole 
states arc calculated to be nearly degenerate in the nu- 
cleus 53 Ca even though there is a noticeable gap in the 
calculated single-particle spectrum (c.f., Fig. [3]). 

To analyze the influence of nuclear deformation on 
the shell evolution, we evaluate the potential energy sur- 
faces of the nucleus 32 Mg which has been intensively dis- 
cussed recently [3814421 ] . Relativistic and non-relativisitic 
mean-field as well as Woods-Saxon calculations result in a 
spherical shape for the ground state of this nucleus [-431 - 
|45| . This problem may be related to the fact that the 
predicted TV = 20 gap is rather large [46| . Whereas the 
observed large B(E2) value for 32 Mg indicates that the 
nucleus is strongly deformed with 02 ~ 0.5 (47[. Our cal- 



FIG. 8. (Color online) PES for the ground states in nuclei 
34 Si (left) and 42 Si (right). 



culations with the re-fitted parameter are plotted in Fig. 
[7] The minimum is around 02 = 0.38 and 7 = 4.8° for 
calculations with kso = ~ K - It should also be mentioned 
that the shape of 32 Mg is rather soft against 02 deforma- 
tion. The nucleus is calculated to be a rigid sphere if we 
take Kso = K - Calcualtions with other parameters [l4j 
lead to a similar conclusion. 

In a recent paper the second + state in 34 Si was ob- 
served, which shows a large deformation parameter of 
02 = 0.29 [H| . The ground state of this nucleus is 
calculated to be spherical with a coexisting shallow de- 
formed minimum (or more exactly a shoulder). The cal- 
culated deformation of the second minimum is around 
02 = 0.38. This second minimum in 34 Si disappears if 
we took kso = K The nucleus 42 Si is calculated to be of 
oblate shape. This agrees with the shell model calcula- 
tions with tensor force shown in Ref. J8j and the rela- 
tivistic mean- field calculations in Ref. [491 ] . The ground 
state in the TV = 28 isotone 44 Si is calculated to be of 
oblate shape with 02 = 0.27 but the minimum is much 
shallower than that in 42 Si. It is expected that in this nu- 
cleus both deformed and spherical congurations coexist 
and mix weakly with each other [Toj | . 

It is expected that the quadrupole collectivity increases 
in neutron-rich nuclei around TV = 40 [50l - [53l |. This is 
supported by our calculations with kso = ~ K i as can 
be seen in Fig. [5] The calculated nuclear deformation 
in 64 Cr is close to that given by the five-dimensional 
quadrupole collective Hamiltonian calculation in Ref. 
[531 ]. Such an enhanced collectivity also indicate that 
TV = 40 do not emerge as a magic number in neutron- 
rich nuclei. 

In summary, we analyzed the shell structure of 
neutron-rich nuclei from a simple phcnomenological 
single-particle point of view. If the SO splitting is rela- 
tively enhanced (i.e., with a strong positive value of kso) 
in neutron-rich nuclei, both HO and Woods-Saxon calcu- 
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FIG. 9. (Color online) PES for N = 40 isotones 64 Cr (left) 
and 66 Fe (right). These nuclei are both calculated to be rigid 
spheres if we took kso = fil- 
iations show that it will destroy the HO magic numbers 
N = 8 and 20 and generate new SO magic numbers like 
TV = 6, 14, 16, 32 and 34. The traditional SO magic num- 
bers N = 28 and 50 and the magic number N = 14 will 



be eroded somehow in neutron-rich nuclei due to the sen- 
sitivity of largcr-Z orbitals to the central potential depth. 
These SO shell closures are more robust than the HO 
magic numbers since their erosion is retarded by the SO 
splitting. In stable nuclei the ld 5 / 2 an d 0g7/2 an d the 
Ohg/2 and I/7/2 orbitals are close to each other. These 
mechanisms both may split those shells, resulting in new 
shell closures like N — 56 and 90. This is in agreement 
with the Skyrme-Hartree-Fock calculations with tensor 
force shown in Ref. Such a simple scheme may give 
a quick estimation on the bulk properties of the single- 
particle spectra. It may also provide a convenient start- 
ing point for a variety of shell model calculations in the 
continuum (see, e.g., Ref. HH) and to explore the effect 
of the pairing correlation and deformation which may in- 
fluence the shell structure. 
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